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Interfacial properties of colloid-polymer mixtures are examined within an effective one-component 
representation, where the polymer degrees of freedom are traced out, leaving a fluid of colloidal 
particles interacting via polymer-induced depletion forces. Restriction is made to zero, one and two- 
body effective potentials, and a free energy functional is used which treats colloid excluded volume 
correlations within Rosenfeld's Fundamental Measure Theory, and depletion-induced attraction 
within first-order perturbation theory. This functional allows a consistent treatment of both ideal 
and interacting polymers. The theory is applied to surface properties near a hard wall, to the 
depletion interaction between two walls, and to the fluid-fluid interface of demixed colloid-polymer 
mixtures. The results of the present theory compare well with predictions of a fully two-component 
representation of mixtures of colloids and ideal polymers (the Asakura-Oosawa model), and allow a 
systematic investigation of the effects of polymer-polymer interactions on interfacial properties. In 
particular, the wall surface tension is found to be significantly larger for interacting than for ideal 
polymers, while the opposite trend is predicted for the fluid-fluid interfacial tension. 



I. INTRODUCTION 

Mixtures of colloidal particles and non-adsorbing poly- 
mers dispersed in a solvent provide experimentalists and 
theoreticians with a very flexible model system to ex- 
plore the statics and dynamics of phase transitions, in- 
cluding fluid-fluid demixing, crystalization, gelation or 
glass transitions Ij, as well as the interfacial properties 
associated with phase coexistence. The flexibility of this 
model system stems from the fact that its properties can 
be easily tuned by varying, among others, the size ra- 
tio of the two components (e.g. by controlling the de- 
gree of polymerization), their concentrations, or the qual- 
ity of the solvent which determines whether the poly- 
mer coils behave essentially like interacting self-avoiding 
walks (SAW), or more like ideal non-interacting poly- 
mers (under ^-solvent conditions). 

More specifically if, as will be done in this paper, one 
adopts an effective one-component representation (by 
tracing out the polymer degrees of freedom) the resulting 
effective interactions between the colloidal particles, ob- 
tained from the well-known depletion mechanism|2|, are 
easily tuned by varying the above physical parameters: 
the depth of the effective attraction between colloids is es- 
sentially controlled by the polymer concentration, while 
the range depends on the polymer size (radius of gyra- 
tion) . The resulting phase diagrams are very sensitive to 
changes in the depletion-induced pair potential between 
colloids(3|. Recent experimental|J| and theoretical^, @ 
efforts have focused on interfacial and wetting properties 
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of colloid-polymer mixtures, either near solid substrates 
(hard walls) or at fluid-fluid phase coexistence. The den- 
sity profiles near a glass wall and the surface tension at 
fluid-fluid phase coexistence were measured 0] and very 
recently the first direct observation of capillaryfluctua- 
tions at the fluid-fluid interface was reported0- Vari- 
ous versions of density functional theory (DFT) of inho- 
mogeneous fluids have been used to determine the den- 
sity profiles, adsorption and surface tension of colloid- 
polymer mixtures near a hard wall, either in the one- 
component 8] or two-component'^, T^ representations, 
or at the phase-separated fluid-fluid interface pH fl^ . 
Most of the theoretical work so far focused on mixtures of 
hard sphere colloids and ideal (non-interacting) polymers 
within the classic model of Asakura and Oosawa 13J and 
of Vrij [l^ (the AO model) , although some attempts have 
been made to extend the DFT calculations to the case of 
non-ideal (interacting) polymer coils^3E|i pointing to 
very significant differences between the two situations. 

This paper presents a unified DFT description of wall- 
fluid and fluid-fluid interfaces within an effective one- 
component representation. The DFT theory is a pertur- 
bative one with respect to the polymer-induced depletion 
attraction between hard sphere colloids, and applies to 
mixtures of colloids and non-interacting, as well as inter- 
acting polymer coils. Apart from a systematic compari- 
son of the density proflles and surface tensions obtained 
for these two cases, the present theory also leads to the 
flrst estimate of the depletion potential between two walls 
induced by colloid-polymer mixtures. 

The paper is organized as follows. Section ^] briefly 
summarizes the effective one and two-component repre- 
sentations of colloid-polymer mixtures. The DFT formu- 
lation used in this paper is presented in section Hill Re- 
sults for colloid density proflles near a hard wall in the 
presence of ideal or interacting polymers are presented in 
section Hvl while the resulting wall surface tension of the 
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mixtures is calculated in section^ The depletion poten- 
tial between hard walls, induced by colloid-polymer mix- 
tures is described in section lVll The fluid-fluid interfacial 
properties calculated within the same DFT approxima- 
tion are presented in section IVIIL while conclusions are 
drawn in section IVIIII 



II. EFFECTIVE ONE AND TWO-COMPONENT 
REPRESENTATIONS 

Consider a binary mixture of Nc colloidal particles and 
Np polymer coils in an external field (e.g. the confining 
field of a hard wall) . The total Hamiltonian of the system 
may be written as 



He 



(1) 



where the colloid-colloid, polymer-polymer and colloid- 
polymer terms JJcc, Hpp and H^p will be assumed to be 
sums of pair interactions between the centres of mass 
(CM) {Ri} {I < i < Nc) and {fj} (1 < j < Np) of 
the colloids and polymers, respectively. Within this as- 
sumption, the individual monomer degrees of freedom of 
the polymer coils have been traced out, and the resulting 
effective pair potential between polymer CM's is a state- 
dependent free energy [T^: similarly the colloid-polymer 
pair potential is a state-dependent free energy resulting 
from a statistical average over monomer degrees of free- 
dom for a fixed distance between the CM's of the two 
particles^. Hence, 



Hcc — Vcc{Rij) 



1<J 

Vpp {Vij ) 

i<3 



(2) 



Hep = J2J2'"'=p^\^' 



while 



Na 



Np 



(3) 



where ip^ and tpp are direct external interactions acting 
on colloidal particles and polymers, respectively. Note 
that (Pp{r) is once more an effective potential acting on 
the polymer CM. 

The colloid-colloid pair potential will henceforth al- 
ways be considered to be of the hard sphere form. The 
AO model for non-interacting polymers further assumes 
that 



Vaflir) 



OO r < (Taf} 
r > (Tap 



a, (3 = c or p 



(4) 



with tJcc = 2i?c {Rc being the colloidal radius), <jpp = 
and Gep — Re + Rp, where Rp is the radius of gyration of 
the polymer coils. A generalization of Rosenfeld's Fun- 
damental Measure Theory (FMT) free energy functional 
for additive hard spheres [T^ been worked out for the 
non-additive hard sphere mixture given by eq^ 18], and 
applied to colloid-polymer interfacial properties||5|,|9j. 

An alternative approach is to trace out the polymer 
degrees of freedom to derive an effective Hamiltonian in- 
volving only the colloid degrees of freedom0,0| . Within 
the semi-grand canonical ensemble, with fixed Ne and 
fixed chemical potential jip of the polymers (or equiv- 
alently, fixed number density Pp of the polymers in a 
reservoir), the effective Hamiltonian is: 



He. 



(5) 



where f2 is the grand potential of the inhomogeneous fluid 
of polymers, which depends parametrically on the posi- 
tions {Ri} of the colloidal particles, can be systemati- 
cally be expanded in terms corresponding to the number 
of colloids HQ: 



Na 



(6) 



n=0 



f2o is the free energy of the pure polymer solution at 
chemical potential ^p, in the volume V of interest, f2i is 
the free energy cost of inserting independent single col- 
loids, takes pairwise colloid-colloid interactions into 
account, and so forth for higher order terms. Such an 
expansion can also be carried out near a single wall of 
area A, where the first few terms are: 



^0 = -Pp{f^p)y +iw,p{f^p)A 



Na 



NeiUl{^ip)+Y,^l"i^^ 

i=l 

Na 

^PVeffiRij) 



(7) 



where Lp%ff{z) is the effective wall-colloid depletion po- 
tential induced by the presence of polymers, Pp{fip) is 
the osmotic pressure of the bulk polymer solution, and 
lw,p{fJ'p) is the surface tension induced by the flat wall. 
Surface tensions are defined throughout relative to the 
position of the hard wall, at z = 0. Inserting a single 
colloid into a bulk polymer solution costs a free energy 
^lifJ'p) = /Nc, and Lp'^J ^ ^z) describes the correction 
to that insertion free energy when the colloid is at a dis- 
tance z from the wall. Veff{Rij) is the effective interac- 
tion between two colloids, induced by the polymer solu- 
tion. The next higher order contributions to f2 include a 
a three-body colloid term, a three-body colloid-colloid- 
wall term, etc... These will be ignored in the present 
work. 

Brader et al. systematically worked out this expan- 
sion for an AO mixture near a flat wall, finding Q: 
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prto = -p;{V - RpA) (since pPpifip) = p; and P^^^^ = 
p^Rp for AO particles[23), = Ncp1^tt{Rc + Rpf + 

/3</?Ao(^i), and /3^2 = Y.i<j P^Ao{.Rij)- PfAoizi) 
acts on each colloid independently and takes the form: 

P^Aois) - - ^(2g - ,s)2(3 + q + s)e{s - 2q), (8) 

where rj^ = AnppRp/S is the polymer packing fraction in 
the reservoir, q = Rp/Rc is the size ratio, and 9{x) is 
the Heaviside step function {9{x) = 1, a; < 0; 9{x) = 0, 
a; > 0). VAo{Rij) is the Asakura-Oosawa depletion'^ 
potential between two spheres in a bath of ideal polymers; 
defining now s = {R — 2Rc)/Rc, vao is given by: 



l3vAo{s) = 



16g3 



{2q~sf{lc, + Aq + s)e{s-2q) (9) 



The effective Hamiltonian for the AO model, with Q, 
restricted to the 0, 1 and 2-body terms in eqElis strictly 
valid only for size ratios q < 0.1547; for larger size ra- 
tios, three and more-body effective interactions between 
colloidal particles come into play j8|, but explicit calcu- 
lations of bulk properties show that they do not play a 
significant role for q up to « 1 |^ . Subsequent calcu- 
lations will neglect more than two-body interactions. For 
the AO model, Oq and the bulk part of fii, f^^"''' = N^uji, 
do not affect the interfacial profiles, but do make con- 
tributions to the surface tension 7 as will be shown in 
subsequent sections. 

The AO pair potential (eq|2I) assumes ideal polymers 
to be spherical, so that the width of the depletion layer 
is always equal to 2Rp. However, polymers can "wrap" 
around colloidal spheres, the more so the smaller the ra- 
dius Rc of the latter. Hence, one expects the width of 
the depletion layer around a colloid to shrink as Rc de- 
creases. This effect can be accounted for by following the 
prescription of refs |H,l3i whereby the AO form (eqEJ 
of the effective pair potential holds, but with a renor- 
malized size ratio q' and a renormalized polymer packing 



fraction r/^': 



q = 



rl 
Vp 



1 + ^q + 



1/3 



Vp 



(10) 



For interacting polymers we can also carry out an ex- 
pansion of fl. Accurate expressions for Pp, ^w,p, and 
Vef f (r) are known and have been validated by computer 
simulations for self-avoiding walk (SAW) chains |20ll23| . 
so that rio, 1^1 and directly follow from eq iQ. The 
Derjaguin approximation, which is reasonably accurate 
for an ideal polymer depletant, is expected to remain 
valid in the case of interacting polymers, at least for suffi- 
ciently small q, so that the required wall-colloid effective 
depletion potential ip%^^ {zi) is given by 



where the effective potential between two colloids in- 
duced by interacting polymers is well approximated by 
the semi-empirical form 

Ms) = -^^Ig^ [s - d{pi)] ' e[s - d) (12) 

Here d — DjRc and D is the range of the depiction po- 
tential: 

Because Vmt {t) depends on polymer density only through 
p^p (i.e. ^p), its form does not change when the polymer 
density is inhomogeneous. In contrast to the AO model, 
higher order colloid-colloid and colloid-wall terms are 
relevant for any size ratio. 

The effective one-component representation of colloid- 
polymer mixtures is thus fully defined, both for ideal and 
interacting polymers. The next step is to define the den- 
sity functional appropriate for the description of the in- 
homogeneous effective one-component model. 



III. DENSITY FUNCTIONAL FORMULATION 

Given the effective Hamiltonian specified by eqs[5land 
one may construct an approximate free energy density 
functional to investigate the properties of the inhomoge- 
neous effective one-component system of colloidal par- 
ticles. The latter interacts via a hard sphere repulsions 
(for r < 2Rc) and a polymer-induced depletion attrac- 
tion that will be described within first order perturbation 
theory [2^, which is expected to be accurate for not too 
small values of q (i.e. for sufficiently long-range attrac- 
tions). The intrinsic free-energy functional is then con- 
veniently split into ideal, hard sphere and perturbation 
parts; for the hard core part we adopt the very accurate 
"White Bear" version 2^ of Rosenfeld's "Fundamental 
Measure" (FMT) functional [TtI 

F [pc{r)] = [pc{r)] + Ffmt [Pc{r)] + Fi [p,(f)] (14) 

The ideal contribution is: 

F,d [Pc] - keT J drpc{r) [Hklp{f)) - l] (15) 

where Ac is an irrelevant colloidal length scale. The 
FMT-hard core contribution is of the "weighted density" 
type, namely 



/3<i5mi(s) ~ 2(3Vint{s), 



(11) 



PFpmt[Pc] = y rfrl<i>i(K(f)}) + $2({n,(f)}) 

+ $3({n,(f)})] (16) 
where the rij (r ) are weighted densities of the form 

nj{r) — / u>'^^\r — r ')pc{r ')dr ' (17) 
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The functions and the weight functions w'-^^ (r) are de- 
fined in Appendix A. Finally, the first order perturbation 
term reads: 



Fi[Pc 



(18) 



with 



df'pc{r')v{\r-r '\)ghs{\r-r '\, pc{r,r")) 

(19) 

where v{r) = wao('') or Vintij") is the depletion potential 
induced by ideal (AO) or interacting (int) polymer coils; 
gjjg is the pair distribution function of the homogeneous 
reference hard sphere fiuid, evaluated at an intermediate 
density between the two points f and r ': 



Pu{r)+ Pu{r') 



with 



47ri?3 



Pc{r')dr'. 



(20) 



(21) 



Pu{r) is a smoothed density profile around r; the radius 
Ri, is of the order of i?c, and results are not expected to 
be very sensitive to the precise value of R^. Following 
the earlier experience of ref we have chosen R^, = 
l.QRc- The form 1 191 is the generalisation of the standard 
first order thermodynamic theory |32| | to inhomogeneous 
fiuids. The choices given in eq [201 and [2] have proved 
very adequate in many DFT calculations of fluids near 
hard walls. 

In the homogeneous limit, where the density profile re- 
duces to the bulk density, the free energy given by eg 1141 
goes over to the Helmholtz free energy of the fluid phase 
calculated within first order thermodynamic perturba- 
tion theory 28] , which leads to reasonable phase diagrams 
of colloid-polymer mixtures |22j|. The familiar generalized 
van der Waals mean field approximation amounts to set- 
ting gHs[f) = 9{2Rc — r), which leads of course to a 
considerable simplification of the DFT calculations. 

Keeping in mind that the total effective Hamiltonian of 
the colloids is given by eqs|SHI| with = f7o+f2i+ri2, the 
grand potential functional to be minimized with respect 
to the density profile Pc{r) is the sum of the intrinsic free 
energy functional El and of the contributions from the 
external fields: 



[Pc] = F [pj + / p,(r) [ip,{r) - /^°] df 



bulk 



(22) 



where p'^ is the bulk chemical potential of the colloids, 
while r^o and fi^"""' = NcOJi are defined in eq[7| Since 
the latter are constants, they have no infiuence on the 
equilibrium profile Pc(f), but they will contribute to the 
equilibrium value of the grand potential, and hence to 
the surface tension, (pd'r) = fci^ + f'c^'^i^ is sum 



of the external potential acting directly on the colloidal 
particles ('^^(t^)) and of the polymer-induced depletion 
contribution, Lp'^^f{'r), given by eqlSjor^Jfor ideal and 
interacting polymers respectively. 

The Euler-Lagrange equation of the variational prob- 
lem reads as: 



Sn^jpcjr^] ^ SF[p,{r)] 
Spc{r) Spc{r) 



+ -Pc^O (23) 



F is the sum of ideal (-Fid) and excess {F^j. = FpMT + Fi) 
parts; defining the local excess chemical potential as 



SF,Mr)] SFpMTiPcm , SFMr}] 



5pc{r) 



5pc(r) 



5pc(r) 



(24) 

the Euler-Lagrange equation may be cast in the form: 
fesTln {hlpcir)] + pT{^ = Pc - V>cir) (25) 

so that the density profile satisfies the non-linear equa- 
tion: 

Pe(r) = p° exp{-/3 [ip^f) + pT{r)) - /i^-^] } (26) 

where p^^''^^ is the the excess part of the bulk chemical 
potential of the colloids. The coupled equations 1241 fwith 
FpMT and Fi given by eas [TCI and [TH|) and must be 
solved numerically by standard iterative procedures. 

Since only planar interfaces will be considered, the ex- 
ternal potential, the density profile and local chemical 
potential all depend only on z, the coordinate orthog- 
onal to the plane of the interface. The corresponding 
simplified expression for the local excess chemical poten- 
tial pI'^ [23 is given in Appendix B, together with the 
weight functions w^^-* (z — z') appropriate for the one- 
dimensional problem. 

Once the equilibrium profile Pc{z) has been deter- 
mined, it can be substituted into the expression of the 
functional 1221 to determine the equilibrium grand poten- 
tial Subtraction of the bulk contribution then allows 
the surface tension 7 to be calculated, according to the 
definition: 



7 



n + pv 

A 



(27) 



where P is the total bulk osmotic pressure of the mixture 
and A the total area of the planar interface. 



IV. DENSITY PROFILES NEAR A HARD WALL 

Consider first the case of a colloid-polymer mixture 
near a hard wall placed at z = 0. The density profiles 
were calculated by solving the one-dimensional versions 
of eas 1241 and 1261 (see Appendix B). A systematic com- 
parison is made between profiles for ideal and interacting 
polymers, over a range of colloid and polymer packing 
fractions 77° and rj^ and for size ratios q = 0.34, 0.5, 0.67, 
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FIG. 1: Colloid density as a function of the distance to the 
hard wall for mixtures of colloids and (a) interacting polymers 
and (b) ideal polymers, at different polymer concentrations. 
In both figures the bulk colloid packing fraction and the size 
ratio were set at 7?° = 0.15 and q = 0.67. 



FIG. 2: Effect of the polymer-colloid size ratio, q, on the 
density profile near a hard wall for rj^ — 0.15 and rjp — 0.4 in- 
duced by (a) interacting and (b) ideal polymers. In both pic- 
tures, the density profile of pure one-component hard spheres 
(without depletant) at 77° = 0.15 as been included for com- 
parison. 



0.85 and 1.05. Figuresn]i,b compare the density profiles 
calculated for a fixed colloidal density and several poly- 
mer concentrations, for interacting and ideal polymers. 
As expected, the colloid adsorption at contact increases 
dramatically with increasing polymer concentration, due 
to the enhanced effective depletion-induced attraction to 
the wall. The effect is stronger for ideal polymers, which 
also lead to strong layering at the highest packing fraction 
hinting at a possible layering transition. The enhanced 
adsorption and layering are easily understood, since the 
polymer induced attraction is substantially stronger for 
ideal polymers (cf. eq^ compared to interacting poly- 
mers feas[T^and [TT|) . at the same Rp, q and p'p 20] . The 
structured profile at r]p = 0.4 in Figure QJa is for a ther- 
modynamic state close to the fluid-fiuid binodal (on the 
colloid-poor "vapour side" ) calculated from the same free 
energy model (see ref [l^) giving further support to the 
possibility of a layering transition,30t. .Slj. 



The effect of the range of the depletion attraction is 
illustrated in Figures |3^,b for interacting and ideal poly- 
mers respectively. The trend of the density profiles with 
increasing q is opposite that observed with increasing Tf^. 
The adsorption at contact is now strongest for the small- 
est size ratio. This is essentially due to the fact that 
the wall-colloid depletion potential at contact, (^^-^-^(O), 
drops as the polymer density decreases. Since ~ 
Tfp/q^ , a small increase in the size ratio q implies a big 
reduction in polymer density, and also in Lp'^J^ (fi). Cal- 
culations carried out at higher colloid packing fractions 
show that the effect of polymer-induced wall-colloid at- 
traction on the density profiles is much reduced as rf^ 
increases, due to the predominance of purely excluded 
volume effects under high density conditionsj2^ [s^ 
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V. WALL SURFACE TENSION 

Once the colloid density profiles have been calculated, 
the surface tension may be determined by substituting 
Pc{z) into eqE21 and then applying relation |^ Some 
care has to be taken in treating the bulk zero and one- 
body contribution SIq and fij""^. These also contribute to 
the bulk properties (calculated by letting Pc{z) go to the 
bulk value p^), and in particular to the colloid chemical 
potential /i^ and to the osmotic pressure of the mixture. 
Thus, for the chemical potential of the bulk one obtains: 



(28) 



where uji is already defined in eq|7| 

In a similar fashion, the total osmotic pressure is the 
sum of the colloid contribution Pc{p^,Pp) (calculated by 
differentiating the homogeneous limit of the intrinsic free 
energy [T^ with respect to p^) and of the osmotic pressure 
of the polymer reservoir: 



P = Pe(Pc,pP+Pp(pP 

Hence, 

n + pv no + ni + n2 + PV 



(29) 



A 

= { 



-PpV + j^^pA 
A 



A 
} + 



(30) 



+ {uJiiPp) J dzpc{z)+ I dzpc{z)(pc{z)} + 
-F[pciz)] n 



+ {- 



A 



/PV 
dzpciz)} + — 



where F [pc{z)] is the colloid contribution to the inhomo- 
geneous free energy, given by eq^^ Substituting eqsEHl 
and 1291 shows that the bulk contributions to flo and fli 
are exactly cancelled so that the final expression for the 
surface tension is given by: 



7^ = lim 

L — 'oc 



lim 

L — ^oo 



n + pv 

A 



(31) 



dz{f [pciz)] + Pc{z)[ipe{z) ~ p,{pl)\ + P,} 



where / — F /V is the intrinsic colloid free energy density 
1141 In other words, the surface tension may be written 
as the sum of a colloid contribution ^w,c (which would be 
obtained by ignoring fio and $7^"""' in eg I22|l and of the 
contribution of the polymers under reservoir conditions: 



(32) 



Equation 1321 is exact within an effective one- 
component representation of a mixture of colloids and 
polymers where in eq |51 is truncated at For 
lw,c{PcT Pp) of interacting (SAW) polymers, we use the 
results of ref . The resulting total surface tensions for 
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FIG. 3: Wall surface tension for ideal (squares) and interact- 
ing polymers (triangles) in the one-component representation, 
and for ideal polymers in the two component representation 
(solid line), as a function of rjp {rfl = 0.15 and q = 0.67). The 
dotted horizontal line is for hard sphere colloids only (without 
depletant). 



the cases of ideal and interacting polymers are plotted 
as a function of the polymer reservoir packing fraction 
in Figure 131 and compared to the predictions of the two- 
component representation (with ideal polymers), calcu- 
lated within scaled particle theorv[lol|. modified so as to 
account for the "polymer wrapping" effect, according to 
the prescriptions of Eq. ^| 



The agreement between the surface tensions calculated 
within the effective one and two-component representa- 
tions is seen to be good for low densities and to deterio- 
rate somewhat beyond. The surface tension with inter- 
acting polymers is substantially larger (up to a factor of 
two at the highest polymer concentrations) compared to 
the ideal polymer case. This is because the colloid-wall 
surface tension is higher with interacting polymers, due 
to the weaker effective wall-colloid attraction, and the 
wall-polymer surface tension is also higher for interact- 
ing polvmerspoj. 

Figure 01 illustrates the variation of the surface tension 
with size ratio q. For fixed values of and rfp, the sur- 
face tension is seen to decrease rapidly as q increases, and 
to tend to the result for hard spheres when q > \. Finally 
the variation of the surface tension with colloid packing 
fraction (for fixed q and yyp is shown in Figure [31 The 
agreement between the effective one and two-component 
representations found at low-packing fractions is seen 
to deteriorate rapidly as 77J? increases, thus illustrating 
the break-down of the effective one-component model 
for highly concentrated colloidal suspensions. 



7 



0.5 
0.4 ■ 
0.3 
QC 0.2 
0.1 
0.0 



^ 



4 {2comp.) 

' {1 comp.) 

A r^"'{1comp.) 
—""—r^fil comp.) 
—^—r„J"'{1 comp.) 
;7;=0.15 Tj;=0.2 



0.4 



0.6 



0.8 



1.0 




FIG. 4: Wall surface tension for ideal (squares) and interact- 
ing (triangles) polymers in the one-component representation, 
and for ideal polymers in the two-component representation 
(solid line), as a function of the size ratio, q {rf^ = 0.2 and 
rfp — 0.2). The dotted horizontal line is for pure hard spheres, 
and the colloidal contributions to the wall surface tension, 
7iD,c, are also shown for comparison. 
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FIG. 5: Wall surface tension for ideal (squares) and interact- 
ing (triangles) polymers in the one-component representation, 
and for ideal polymers in the two-component representation 
(solid line), as a function of the bulk colloid packing fraction, 
rf^ {r)l = 0.2 and q = 0.67). 



VI. DEPLETION INTERACTION BETWEEN 
HARD WALLS 

The depletion potential per unit area induced by a 
colloid-polymer mixture between two hard walls of area 
A separated by L is determined by: 



W{L) 



n{L) - n{oo) 



(33) 



Clearly, W[L) ^ as L ^ oo, while W{L = 0) = 
-27m, because when the walls come into contact two 



FIG. 6: Effective depletion potential between two flat parallel 
walls as a function of the wall-wall distance, L, induced by 
pure hard colloids (dotted line), colloids plus ideal polymers 
(triangles) and colloids plus ideal polymers (squares), rf^ = 
0.15, ril = 0.15 and q = 0.67. 



fluid-wall interfaces are destroyed. For arbitrary L, the 
colloid contribution Wc to W{L) obtained by momentar- 
ily ignoring the bulk contributions JIq and ilj"'*^ to the 
total grand potential, reads, by straightforward analogy 
with expression 13 II for the surface tension: 



W,{L) 



dz{F[p,{z)] 



+ pc{z)[ipc{z)- ^ic{pl)]+Pc}-2-i^,c (34) 

Figure shows the comparison of the depletion poten- 
tials Wc{L) calculated for a pure hard sphere depletant 
(i.e. vip — 0) and for hard sphere colloids and ideal or 
interacting polymers, within the present effective one- 
component representation. The extra attraction between 
the walls and colloidal particles forces the latter into the 
region between the two plates, enhancing the potential 
barrier of the depletion potential aX L — R^. The effect is 
stronger for ideal than for interacting polymers, because 
effective wall-colloid and colloid-colloid attractive inter- 
actions are enhanced in going from interacting to ideal 
polymers. For L < 2Rc the hard colloids are excluded 
from the space between the walls, and the depletion po- 
tential is entirely controlled by the external, unbalanced 
osmotic pressure which pushes the walls together, i.e.: 



(35) 



To obtain the total depletion potential, the volume 
terms fto and fi^"""' must be included in the grand po- 
tential. It is easily verified that, within the truncation of 
f2 at r22, i-e. ignoring colloid-colloid-polymer and simi- 
lar higher order terms, the depletion potential takes the 
form: 



W{L) 



Wc{L) 



p;l 



L<Dl 
L>DP 



(36) 
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FIG. 7; Comparison between the wall-wall depletion poten- 
tial in a mixture of colloids and ideal polymers, with and 
without the bulk contributions to the free energy {W and 
Wc, respectively). t]° = 0.15, rip = 0.15 and q — 0.67. 



where Pp — ksTp^ is the osmotic pressure of the poly- 
mers in the reservoir, ^w,p = ksTppRp is the correspond- 
ing wall-polymer surface tension, and = '2"fw.p/Pp 
is the range of the polymer depletion potential, which 
takes the value Z3J^ = 2Rp for AO particles. The result- 



ing total depletion potential for ideal polymers is shown 
in Figure 13 It remains continuous at L 
resulting force is discontinuous at that separation. 



2Rp, but the 



The basic input in the evaluation of surface tensions 
or depletion forces are the density profiles Pc{z). In 
order to check the reliability of the profiles calculated 
within the present perturbation DFT, we have carried 
out some grand canonical Monte-Carlo (GCMC) simula- 
tions of profiles, using the same effective colloid-colloid 
and colloid- wall interactions as with DFT. Two typi- 
cal examples of such a comparison are shown in Fig- 
ure |HK,b for a mixture of colloids and interacting poly- 
mers confined between two hard walls (slit geometry). 
At the lower polymer concentrations, corresponding to 
relatively weak colloid-colloid and wall-colloid attrac- 
tion, the agreement between DFT and simulation is ex- 
cellent (Figure |S^). At the higher polymer concentration 
(Figure inh) the attractive perturbation is much stronger 
and, as expected, the agreement worsens. Nevertheless, 
the perturbation DFT is still capable of reproducing the 
main features of the density profile. FigurelHt also shows 
the result from the mean-field DFT, which amounts to 
neglecting correlations in Fi [pc\ (cf . section 0) . The 
disagreement with the simulated data is now much more 
severe (see the inset in Figure (SId), stressing the impor- 
tance of properly including reference system correlations 
in eqUl 
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FIG. 8: (a) Colloid packing fraction between two hard walls 
obtained using GCMC simulations (squares) and DFT (solid 
line) for a mixture of colloids and interacting polymers and 
for q = 0.67, 77° = 0.15 and rfp = 0.15. (b) Same as FigureHi, 
but with rfp — 0.5. Dashed lines are the results obtained from 
the mean-field DFT theory. The inset shows a magnification 
of the inner part of the data for better comparison. 



VII. THE FLUID-FLUID INTERFACE 

We finally turn to the fluid-fluid interface at coexis- 
tence between colloid-poor ( "gas" ) and colloid-rich ( "liq- 
uid") phases. This is a "free" interface, since transla- 
tional symmetry is broken in the absence of any external 
potential, i.e. V'c(^) = 0. An accurate and consistent 
determination of the phase boundary (binodal) in the 
ryj? (colloid packing fraction)-/^^ (polymer reservoir pack- 
ing fraction) is required. The bulk free energy per unit 
volume, obtained from eqs I14H19I in the homogeneous 
{pc{r) = Pc) limit reads: 



r;0(4-3r?0) 



(1 - 7yO)2 
drgnsir, p°)Pv{r). 



(37) 
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This free energy depends implicitly on the polymer con- 
centration ifp via the effective depletion pair potential 
v{r) = v(r;r]p), which is given by eqEl (vao) for ideal 
polymers and by eg 1121 (vi„-t] for interacting polymers. 

For any fixed value of r]p the colloid packing fractions 
of the coexisting phases may be determined by applying 
the standard Maxwell double-tangent construction to /. 
The resulting binodals for ideal and interacting polymers 
(not shown) are very close to those calculated in ref , 
where an estimate of the second order perturbation cor- 
rection to the free energy was included. The trends of 
the binodals are according to expectation. Qualitatively, 
the behaviour for ideal and interacting polymers is sim- 
ilar, with critical points at higher rjp as q increases, but 
with considerably "flatter" binodals for interacting poly- 
mers. There are, however, very significant quantitative 
differences, since the binodals for interacting polymers 
are shifted to considerably larger values of rj^ for each size 
ratio q, i.e. polymer interactions enhance the miscibility 
of colloid-polymer mixtures. This is easily understood, 
since the polymer-induced depletion attraction between 
colloids is weaker (for given q and ryp for interacting poly- 
mers. 

The colloid density profile at the fluid-ffuid interface 
is calculated by minimizing the grand potential fea I22|l . 
with (fie = 0, which leads back to the exr)ressionl26l f again 
with ifc = 0). The profiles for ideal and interacting poly- 
mers are shown in Figures |H| and 1101 under conditions 
close to the fiuid-fiuid-solid triple point, for a size ra- 
tio q = 0.67. The profiles are compared to earlier pre- 
dictions based on a square gradient density functional 
(SGT) 0, using the same effective one-component de- 
scription. There are remarkable differences between the 
profiles corresponding to interacting and ideal polymers. 
The interface is much sharper in the latter case and ex- 
hibits striking oscillations on the high colloid density 
( "liquid" ) side 9] . The profile obtained with interacting 
polymers varies more smoothly, and shows no sign of os- 
cillations. This clear-cut difference in behaviour may be 
partly understood by noting that the jump in colloid den- 
sity between the two fluid phases is significantly smaller 
in the case of interacting polymers, for which the triple 
and critical points arc much closer than for ideal poly- 
mers. Figures 1^ and [TUl also show the result of our earlier 
square gradient calculations (SGT) The correspond- 
ing profiles agree reasonably well with the results of the 
more elaborated density functional used in this paper; as 
expected the square gradient functional cannot account 
for the oscillatory structuring on the dense fluid side in 
Figure Uni 

In the case of interacting polymers (Figure ^ the 
width of the interface predicted by the present pertur- 
bation DFT is significantly larger than that obtained 
from square gradient theory. The interfacial width w 
is conventionally defined as the distance between the 
two points where the density profile pdz) reaches 90% 
and 10% of the difference between the bulk densities 
of the two coexisting phases. Results for w obtained 



0.40 




FIG. 9; Comparison between the colloid density profile 
for interacting polymers obtained using our perturbative 
DFT (solid line) and the square gradient approximation 
SGT(dashed line), for rjp — 0.83 and q = 0.67 (close to the 
triple point). 




FIG. 10: Comparison between the colloid density profile for 
ideal polymers obtained using our perturbative DFT (solid 
line) and the square gradient approximation SGT (dashed 
line), for rjp = 0.56 and q — 0.67 (close to the triple point). 

within the perturbation DFT are plotted in Figure ITTI 
as a function of the deviation of the polymer reser- 
voir packing fraction from its value at the critical point, 
^ — iVp ~ Vp'^^^*) I Vp'^"^ ^ for ideal and interacting poly- 
mers and for three size ratios q. As expected, the width 
increases with q and increases sharply as a decreases, i.e. 
upon approaching the critical point. Although the widths 
w for ideal and interacting polymers are rather close for 
given a and g, polymer interactions tend to reduce the 
width compared to the ideal found in Il3|- 

The surface tension 7 of the fluid-fluid interface is 
given by a relation similar to eg 1311 and it is easily ver- 
ified that the bulk terms fio and rj^"'*^ in eq|221do not 
contribute to 7. The relation reads: 

/oo 
dz{f[p,{z)]- p,{z)^il + P,} (38) 
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FIG. 11: Interfacial width obtained using the perturbation 
DFT for interacting polymers (solid lines and black symbols) 
and ideal polymers (dashed lines and white symbols), and for 
q = 0.34, 0.67 and 1.05. 
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FIG. 13: Surface tension of the free liquid-gas interface for 
interacting polymers (solid lines and black symbols) and ideal 
polymers (dashed lines and white symbols), and for q — 0.34, 
0.67 and 1.05. 
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FIG. 12: Comparison between the surface tension of the gas- 
liquid interface obtained from our perturbation DFT (solid 
lines and black symbols) and the square gradient theory SGT 
(dashed lines and white symbols), for ideal polymers. 



mers is finally made in Figure 1131 The surface tension 
is seen to be lowered when polymer interactions are in- 
cluded, a trend which is opposite to that observed for 
the wall surface tension -fw (cf. Figure 01 . A fluid- fluid 
interface also implies an inhomogeneous polymer density 
profile. For ideal polymers at constant fip this does not 
cost any additional free-energy, but this is not true for in- 
teracting polymers. For the latter depletant, this correc- 
tion will raise the surface tension compared to the values 
we calculate within the one-component picture. On the 
other hand, for ideal polymers, the many-body interac- 
tions between the colloids, fully taken into account with 
the two-component formulation, also result in a larger 
surface tension. The effect of many-body interactions for 
interacting polymers is unknown. Thus our prediction 
for the relative strength of the surface tensions should be 
tempered by the fact that the exact form of the neglected 
corrections is not know and awaits a full two-component 
treatment of interacting polymer-colloid mixtures. 



where /ij; and Pc are the common values of the colloid 
contribution to the chemical potential and of the pres- 
sure in the coexisting bulk phases. Figure [T^ compares 
results obtained from the perturbation and square gradi- 
ent theories for ideal polymers, at three different size ra- 
tios; 7 is plotted versus the difference between the colloid 
packing fractions in the "liquid" and "gas" phases. The 
surface tensions are seen to increase with q and tj^ — tj^ ■ 
The agreement between the two theories is excellent for 
q = 1.05, and deteriorates for smaller size ratios. The 
surface tensions obtained from perturbation DFT are 
systematically larger than their square gradient counter- 
parts and closer to the results of the two-component DFT 
of ref 0. 

A comparison between the surface tensions calculated 
within perturbation DFT for ideal and interacting poly- 



VIII. CONCLUSION 

We have introduced a simple generic density functional 
description of colloid-polymer mixtures within the effec- 
tive one-component picture where polymer degrees of 
freedom have been traced out. The free energy func- 
tional accounts correctly for the colloid excluded volume 
effects, and treats the polymer-induced depletion attrac- 
tion within first order perturbation theory. The DFT for- 
mulation is very fiexible and can be applied to confined 
colloid-polymer mixtures as well as to the free interface 
between coexisting fluid phases. Although this effective 
one-component description accounts only for pair-wise 
additive depletion interactions, and neglects more-than- 
two-body effective interactions (which are automatically 
included in an effective two-component representation). 
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the present theory has the advantage of being able to 
treat the cases of ideal as well as interacting polymers 
consistently, thus allowing a direct estimate of the ef- 
fect of polymer interactions on the interfacial properties 
of colloid-polymer mixtures. The present effective one- 
component DFT description can be tested against the 
more fundamental two-component DFT picture in the 
case of ideal polymers^ J^. Reasonable agreement is 
found for size ratios up to g « 1, giving confidence in the 
predictions of the simpler one-component representation 
for ideal and interacting polymers alike. The reduction 
to an effective one-component representation leads to the 
appearance of zero and one-body bulk contributions to 
the total grand potential of the mixture (terms Qq and 
flbuik jj-^ eg I22|l . which make significant contributions to 
the wall surface tension 7^, and to the depletion interac- 
tion between two walls, and must not be overlooked. The 
key findings of the present investigation may be summa- 
rized as follows: 

• Density profiles Pc{z) near a hard wall deviate con- 
siderably from pure hard sphere behaviour at low 
colloid packing fractions. Significant differences be- 
tween the behaviour observed for ideal and inter- 
acting polymers may be traced back to the weaker 
depletion attraction induced by the latter. 

• When the bulk terms flo and fij""^ in the grand 
potential are properly taken into account, the wall 
surface tensions 7^, calculated with ideal polymers 
are in good agreement with a recent full two- 
component treatment of the AO model 10], for 
g < 1 and not too high colloid packing fraction 

Vc- 

• The depletion potential induced by colloid-polymer 
mixtures between two walls exhibits structure on 
both length scales Rc and Rp which opens the pos- 
sibility of further flexibility in "engineering" effec- 
tive depletion forces. 

• In general, the wall surface tension induced by 
mixtures of colloids and interacting polymers is 
larger than that found for ideal polymers. This 
trend is opposite to that predicted for the fluid- 
fluid interfacial tension. The present predictions 
agree at least qualitatively with the experimental 
data of ref 0, and may be understood in terms of 
the wall surface tension of the pure polymer com- 
ponent. 

• The present effective one-component DFT predicts 
an oscillatory density profile on the high colloid 
density side of the fluid-fluid interface in the case of 
ideal polymers sufficiently close to the triple point; 
this behaviour agrees with the earlier prediction 
based on the effective two-component (AO model) 
representation The oscillations are not observed 
in the case of interacting polymers, presumably due 



to the smoother variation (larger width) of the in- 
terfacial profile. 

• A direct comparison between the results of square 
gradient analysis^^ the present more sophis- 
ticated DFT shows that the former is surprisingly 
accurate, even for rather sharp fluid-fiuid interfaces 
(except, of course, as regards the oscillatory be- 
haviour observed for ideal polymers). The good 
agreement validates the predictions of the sim- 
ple square gradient theory concerning the effect of 
polymer-polymer interactions on the fluid-fluid in- 
terfacial properties^]. 

The present effective one-component description pre- 
dicts a number of signiflcant differences between inter- 
facial behaviour of mixtures of colloids and ideal ver- 
sus interacting polymer coils. The corresponding two- 
component representation involving ideal polymers (i.e. 
the AO model) is well established, and conflrms many 
of the predictions of the present effective one-component 
picture. Future work should focus on developing a viable 
effective two-component representation in the case of in- 
teracting polymers, in order to validate the present pre- 
dictions based on the effective one-component descrip- 
tion. 
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IX. APPENDIX 

A. Fundamental Measure Theory for hard spheres. 

In order to model the hard-sphere nature of the col- 
loidal particles, we have used the White-Bear version of 
the Rosenfeld functionalfl^ . This functional improves 
the one proposed by Rosenfeld^3] in that it leads to 
the Mansoori et al. equation of state for a homoge- 
neous mixture [3^. The free energy density ^fmt — 
$1 + $2 + *3 readslii]: 



$1 - 
$2 = 



-no ln(l - ns) 

1-^3 
(nl - 3n2ny2, 



(39) 



- (1 - ?i3)^ln(l - 713) 
367rn§(l - ns)^ 
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where {nj{r)} are weighted densities, obtained as convo- where v{r) is the corresponding colloid-colloid depletion 
lutions of the colloid density pdr) and the weight func- potential, its range and 
tions w*^^-* (see ea I17II . The latter are given by: 



UJ^^^r) = e{Rc-r) 

(f) = (f/r)(5(i?c - r) 



,iV2) 



47ri?2 



47ri?^ 



(f) 



(40) 
u;(^2)(f) 
47ri?^ 



^I^^=c.(^)(z'-z) 



(42) 



B. Particular case: planar geometry. 



For planar geometry, the weight functions are given by: 



For the case of planar fluid-fluid interfaces or infi- 
nite planar walls, the calculations simplify to a one- 
dimensional problem. Then, the position-dependent ex- 
cess chemical potential ij,'^^{r), which can be obtained 
from the functional derivative of the excess part of our 
perturbation free energy functional (ea l24|l . depends only 
on z, the distance to the wall (or to the interface). Its 
explicit form for z — zi ispfll]: 



PpT(z,) 



■\p d^FMT [njiz2)] ^^j^^'^) 
^ dnj{z2) Spc{zi) ^ 



(41) 



27r 



dz2Pc{z2) X 



+ 



\zi~Z2\ 

37r 



dri2ri2l3v{ri2)ghs{ri2, Pc{zi, Z2)) 



dz2Pc{z2) [RI - {Zl - Z2f] X 



dzspciz^) X 



\Z2-Z3\ 



dr23r23f3v{ri2) 



dghsir23,Pciz2,Z3)) 



J'^z' 


— z) 






J'\z' 


— z) 


= U{R,~\z-z'\) 




J'\z' 


— z) 

— z) 


= 2TrRc0{Rc~ \z-z 'D 
= i:[Rl-{z-z 'f] 6{R, 


— \z — z 


uj^^'\z' 


— z) 

— z) 


= ^^{z-z')9{R,,~\z- 
= 2tt{z-z ')e{R^- \z-z 


z'\)k 

'l)fc 



The ^-dependent colloid density profile is obtained as 



Pciz) = p°exp{-/3 [Mz) + pTiz)) - M^'^ll (43) 



The coupled equations eas 1411 and 1431 are solved itera- 
tively. 
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